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PHYSlC/iL AND HATIILMATIt'AL MODKLLING OF LADLE METALLURGY OPERATIONS 

N. El-Kaddah and J. Szekply 
Massachusetts Inst’tute of Tochnoloqy 
Cambridge, Massachusetts 02139 USA 

ABSTRACT 

Experimental measurements are rei>orted, on the velocity fields and 
turbulence- parameters on a water model of an argon stirred ladle. These 
velocity measurements are complemented by direct heat transfer measure- 
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ments, obtained by studying the rate at which icc rods- immersed into tlie 

ei 

\ system melt, at various locations. The theoretical work undertaken involv- 

n 

B 

cd the use of the turbulence Wavier-Stokes equations in conjunction with the 
k-E model to predict the local velocity fields and the maps of the turbu- 
lence parameters. These theoretical predictions were in reasonably good 
agreement with the experimentally measured velocity fields; the agreement 
betwi.'on the predicted and the measured turbulence parameters was less per- 
fect, but still satisfactory. The implications of those findings to the 
modelling of ladle metallurgical operations arc discussed. 

1 . INTRODUCTION 

During the past decade there has been a growing interest in the develop- 
ment and application of ladle metallurgy as a means of the finlsiiing treat- 
ment of molten steel. Indeed some form of a ladle metallurgy step, such as 

argon stirring in order to achieve homogenization has become an established 

(1 ) 

■ component of the- overall stoclmaking sequence, ' 

Parallel to the industrial develoments a great deal of useful research 
has been done on various aspects of ladle metallurgy. Ic is generally 
appreciated that apart from process chemistry the fl-iid flow phenomena, such 
.as the flow patterns and the extent of agitation play a key role in deter- 
mining tlm overall cffic ->/ those opt'srations. In gen*iral ladle 
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metallurqical systems may be agitated cither by the injection of gas 

streams or by the use of an electromagnetic force field. In view of the 

topic of this conference our attention will be restricted to gas bubble 

driven circulation systems only, but within this framework we shall seek 

to review recent efforts aimed at the better understanding of the behavior 

of these systems through the use of 

mathematical modelling 

physical modelling 

plant scale experimentation. 

In considering gas bubble driven circulation systems in general, such 
as an argon stirred ladle, from the standixsint of transport phenomena we 
are faced with three groups of problems: 

(i) The definition of the gross features of the circulation system, 

that is the overall flow patterns and the mixing time. ^ 

(ii) The definition of heat mass and momentum transfer between the 
melt and the bounding surfaces, i.e. problems concerned v;ith refractory 
erosion, slag - metal reactions and oxygen transfer from the refractory 
to the bath and finally 

(iii) The description of the detailed structure of the turbulent flow, 
such as the nap of the turbulent energy dissipation which would be needed 
to model deoxidation kinetics, the coalescence of inclusion particles 
. and in general the interaction between suspended particles and the melt. 

In reviewing previous work, major advances have been made regarding 
the definition of the overall flow patterns and the mixing times, through 

the use of physical modelling and the application of various turbulence 
(*— 13 ) 

models. Indeed it may be stated with some confidence that these problems 

are well uncersttiod for certain geometries would be essentially a routine 


task 
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• As an illustration Piq. 1 shows a plot of the mixing time against i 
the rate of energy dissipation (i.c. the total ix>wcr input into the 
system) and it is seen that most systems of metallurgical interest may 
be reasonably well represented on this plot. 

It should be noted that tticre are sound fundamental reasons why 
such an improved understanding of these systems was readily developed. 

As far as the circulation of these systems is concerned in the bulk the 
inertial terms represent the key factors, thus the overall flow field 
may be readily modelled (either mathematically or physically) . /en if 
the viscous (including the eddy phenomena) components are not accurately 
represented. This is a key feature of body force driven flows. 

The situation is rather loss satisfactory regarding the other two 
problem areas, namely fluid to wall (or slag - metal) trans^)ort and the 
definition of the fine structure of the flow. Notwithstanding their 
obvious practical imiiortanco, these aspects of the problem received very 
much less attention; while attempts have been made to extend the theo- 
retical treatment to the explicit representation of these phenomena this 
work has to be regarded as largely tentative at this stage, because of 
the absence of hard corroborating evidence. 

The purpose of the work to be described in this paper is to report 
on recent research addressed to these problems. The actual results will 
be presented under the following tliree headings: 

(1) Mathematical Model Pcveloi’ment 

(2) Physical Modelling and tlic Comparison of the 
Measurcment.s with Predictions 


•*A 
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(3) Plant Scale .Measurements 
2. MATHEMATICAL MODELLING 

In essence the qas bubble driven circulation system represents 
an axi-symmetrical turbulent recirculating flow problem of the type 
which has been tackled by numerous investigators, One 

of the main problems encountered in previous studies was to obtain 
a proper representation of the boundary between the qas bubble rich 
jet cone region and the bulk of the liquid, both in*tqrms of 
position and in terms of stating the proper boundary conditions for 
the velocity and the velocity gradients. 

The following assumptions have been made in the developing 
of the fluid flow equations: 

(1) The two phase region was confined to a plume, the 
dimensions of which were determined experimentally, e.g. from the visual 
observation of the "break through zone". 

(2) The two phase region was assumed to be homogeneous 

(3) The velocity and the momentum flux were assumed to be contin- 
uous across the jet cone (plume) boundary 

(4) The two dimensional, turbulent Navier-Stokes equations were 
used in conjunction with the k-c model to represent the fluid flow field, 
w)iich was considered to bo driven by the density difference in the two 
phase region. 

Assumptions (l)-(3) arc illustrated in Fig. (2), 

Within the framework of these assumptions the governing equations take 


the following form: 






equation of contunity: 


r I? ■ “ 


( 1 ) 


Momentum balance in the z-direction 
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( 2 ) 


Momentum balance in the r-direction 
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where 


(3) 


P = , r > 


(4) 


p = a. p + (l-a)p,, r<r 
g . x. c 


(5) 


It follows that for the purpose of modelling the whole domain 

has been treated as a homogenous medium, but with a spatially variable 

density. Using drift flux model which allows the slip J>otwoon the 

(14) 

bubbles and the fluid, the quantity ct may be expressed as 


2TI 


Q - vr^ a(l-a)U 
g c " 

ru dr 
o z 


( 6 ) 
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where is the rising velocity of a characteristic single bubble, which 
is estimated as 40 cm/s. 

The boundary conditions for the equation of motion take the 
following form: 
at the axis 



at the walls 


(7) 

( 8 ) 


» • 

U = U = 0 . . (9) 

r z 

at the free surface 


0, = 0 ; z=H (10) 

z 

du ' 

^ = 0 ; z = H (11) 

at the orifice 

U=U;r = 0,z = 0 
z o 

(12) 

or U=0 ; r = 0 , z>0 

z 

Eqs. (1-12) constitute the complete statement of the fluid flow 

problem except for the definition of effective viscosity. In the present 

case we used the well cstablis)icd k-C model for computing the turbulent 
(15,16) 


viscosity. 
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3. PHYSICAL MODELLING AND THE COMPARISON OF MKASUREMKNTS WITH 
PREDICTIONS . 

3.1 The apparatus 

The apparatus was so constructed to allow the convenient measure- 
ment of the velocity fields and the turbulence characteristics in a 
cylindrical tank containing water, which is being agifated by an 
ascending gas bubble stream. 

• . 

In essence the apparatus consisted of .a cylindrical tank, contain- 
ing water which was agitated by a gas stream, introduced through the 
bottom of the container, via a centrally located orifice. This cylin- 
drical tank was surrounded by a square tank, also filled with water, 
in order to eliminate the parallaxis effects in the optical measure- 
ments. 

Fig. 3 shows a schematic sketch of the apparatus, which also indicates 
the principal dimensions of the container. Fig. 3 also shows the support- 
ing and positioning assembly used for the melting of ice rods. The actual 
melting rate was determined using photographic technique. 

A schematic outline of the apparatus as a whole is shown in Fig. 4, 
where it is seen that the measurements of the velocities and of the 
turbulence characteristics was accomplished through tlie use of a laser 
anemometer. A TSI Laser nnomometer was used, employing a Spectrophysics 
No. 124, lle-Nc laser. The anemometer was operated in a dual beeun mode. 
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with a analysis of the forward scattered light. 

The velocity measurements involved the use of essentially standard 
procedures commonly employed in laser anomometry. The Reynolds stresses 
were obtained by taking measurements regarding a certain direction and 
then displacing these by 45° following the relation: 

^ = \ - U'^l . (15) 

2 2 

where Uj and U' ..u ^ • .. .^cO ^ ._o , 

1 2 are the fluctuating velocities at +45 and -45 from 

the coordinates i and j. 

The procedure for obtaining Reynolds stresses by such means is 
again well documented. 

3.2 Experimental Measurements and Computed Results 
In this section we shall present a selection of the experimental 
measurements and these will, be compared with the computed values of the 
corresponding parameters. The material to be presented will include 
information on the velocities, on the turbulent kinetic energy and on the 
Reynolds stresses. 

The mean Velocity Field 

Figs. 5a, b show a comparison between the experimentally measured and 
the theoretically predicted velocity fields. It is seen that the agree- 
ment between moasuremonts and predictions is reasonably good. It should be 
noted that the use of k-C model in conjunction with approi>riate wall 
functions will give a significantly better agreement between the measure- 


ments and the predictions, especially in the near wall regions than would 
the use of simpler computational procedures, employing a single value 
of the effective viscosity. 
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Figs. 6a,b,c shows a comparison between pf<>dictcd and the measured 
absolute values of the velocity vector at three axial {positions. It is 
seen that the agreement is quite reasonable throughout the domain. 

rig. 7 shows the computed velocity distribution in the jet cone 

region, which has been normalized with respect to the maximum (centerline) 

velocity. It is seen that this normalized axial velocity is essentially 

independent of the axial position. The Gaussian type distribution 

found seems to be in good agreement with the experimental measurements 

(9) 

of Tse-Chang et al. It is noted that measurements were also made of 

Ug which we found to be close to zero in all cases. However, as will 

he shown subsequently Ul* f 0, 

t) 

TURBULENCE CtlARACTKRISTICS OF THE SYSTEM 

Fig. 8 shows the profiles of the rros velocity components normalized 
with respect to the mean velocity in the jet cone at the corresponding 
vertical position. A somewhat simplified analytical method was used 
for estimating the latter quantity. The turbulence appears to be 

fairly isotropic except in the vicinity of the solid surfaces (side 
walls andthe bottom), which is essentially in line with expectations. 

The other interesting feature of these results is the similarity exhibited 
by the plots shown in Figs. 9a, b. 

Figs. 10a, c show a compari.son between the cpxerimentally measured 
and the calculated values of the turbulent kinetic energy. The agreement 
seems to be within about a factor of two in most cases, which is not 
unreasonable; but this discrepancy indicates that further study of these 
phenomena would bo fully justified. This further theoretical work could 
involve .’ither the adjustment of the parameters in the k-e model or perhaps 
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3.3 HEAT TRANSFER AND PI I ASK CHANGE 

In these series of experiments we studied the rate at which ice 
rods melt, upoh immersion into a gas bubble agitated water |x>ol. The 
main motivation for this work was to obtain a better insight into the 
effect of the velocity fields and the local turbulence parameters in 
the system on heat transfer between the fluid and an immersed solid 
body. These types of problems are of course of considerable practical 
interes*^ both regai'ding wall erosion and dissolution kinetic. 

Fi^. 12 show typical photographs of the molting ice rods for 
various values of time. , , 

Fig. 13, show typical plots of the position of the melting inter- 
face against time, obtained within the jet cone and outside the jet 
cone respectively. The near linearity of those plots indicates that the 
time required for the attainment of steady state melting conditions was 
relatively short. Figs. 14-18 show plots of the initial melting rate as a 
function of position within the agitated liquid. 

Figs. 14 and 15 show the radial variation in the melting rate of 
two constant axial positions. It is seen that the melting rate is the 
highest in the jet cone and that it decreases quite rapidly with the 
radial distance from the axis of the gas jut. 

Fig. 16 summarizes the radial variation in the melting rates at 
various axial positions in jet cone. The maximum melting rate is seen to 
occur at the center at an intermediate height; it may be shown that the 
fluid velocity is the maximum at this position., 
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Pigs. 17 and 18 shows the axial variation of the melting rate for two 
fixed radial -positions within the tank, one near the wall and the other 
in an intermediate position. Figs. 17 and 18 show a monotonous increase 
iii the rate of melting with the distance from the bottom of the vessel. A 
sharp increase is observed on approaching the free surface. This behavior 
is attributable to the high level of turbulence in the vicinity of the free 
surface. Pig. 17 depicting the behavior of the system at an intermediate 
radial position, shows a minimum in the melting rate, corresponding to the 
center of the vortex, which of course, is a quiescent portion of the vessel. 

4. PLANT SCALE MEASUREMENTS 

At present there are only a limited number of plant scale measurements 
where sufficient details are available on the system to allow a meaningful 
comparison between the measurements and the predictions. Nonetheless it is 
important to cite these because they provide the ultimate test of the results 
obtained in the mathematical and water modelling work wltich has been described 
in the bulk of the paper. 

4.1 Mixing 

Mixing is the topic which has been the most extensively studied. Indeed 
numerous data are available on tracer dispersion and mixing a vessel of 
various size. Since this field is now quite well understood, wo 

shall present only one, typical example, illustrating that the theoretical 
predictions provide quite satisfactory answers regarding the mixing time. It 
should be noted that in most industrial scale systems, where the flow is 
reasonable turbulent, the dispersion of the tracer takes place by two 
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mechanisms, namely bulk convection and eddy diffusion. The actual 

convective time scale is very short, in these systems, typically of the 

order of a few seconds (linear dimensions being of the order of 1-2 m, 

while the characteristic velocity being about 0.3 m/s). Thus the rate 

determinine step in tracer dispersion or mixing will bo eddy diffusion 

between the respective streamlines with eddies of smaller time scale 

1/2 

(Kolmogoroff eddies t=(vi/c) ~*ls) as sketched in' Fiq. 19. This behavior 
is of course consistent with the discussion presented in the introductory 
section, concerning the fact that the mixing characteristics of a system 
may be readily related to the total energy input. 

4.2 Desulfurization Kinetics 

Fig. 20 shows a sketch of a typical arrangement used for the desul- 

( 18) 

furization of steel melts. It is seen that a soli i desulfurization 

agent is injected_ into the melt together with a carrier gas, through a 
suitably placed lance. The overall rate of the destJlfur ization reactions 
depends on two factors; in addition to the obvious criteria determined by 
thermodynamic equilibria. 

(i) The rate at which the individual particles react is determined 

(19 21) 

by the Kilmogoroff eddies, ' which clearly depends on the local rate 
of turbulent energy dissipation. 

(ii) The rate at which the reactant, that is the sulfur is being 
transiiortcd to the regions where the reaction is actually taking place. 

hs mentioned earlier, in connection with mixing, tlic rate of this transport 
will be limited by eddy diffusion. 
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Figs. 21a, b, shows a comparison between the experimentally measured 

and the theoretically predicted desulfurization rates, for a pilot scale 

(18) 

and an industrial scale system respectively. It should be remarked, 

that while these two sets of curves appear to be juite similar the actual 
reaction mechanism differed substantially in the pilot scale and in the 
full scale system. 

In the pilot scale operation, wiiere the system was heavily deoxidized, 

t 

the reaction took place in the plume, while in the industrial scale unit 

which was not strongly deoxidized, most of the desulfurization took place 

• > 

at the slag metal interface. The very good .mixing which was provided by 
the gas bubbles in the large scale system provided an adequate rate of 
transport to the reactions zone. 

4.3 Moat or mass Transfer to Imn'orsed Surfaces 
A key problem in the opiation of ladle metallurgical operations 
is the definition of the factors that govern wall erosion, and hence the 
lining. Relatively little modelling work has been done in thin area, 
but a useful start has been made regarding the study of the rate at which 
immers.:d bodies dissolve. 

Fig. 22 shows a typical experimental arrangement, used at MEFOS in 
the study of the rate at which immersed graphite rods would dissolve in 
6 ton iron melts, and the actual experimental results are given in Fig. 23 
in the form of the local mass transfer coefficients. 

The numerical values of these mass transfer coefficients would be 
difficult to predict from first principles (and the same applies to the wall 
erosion rates) because the local fluxes at the solid surface will depend 
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not only on the local velocities iHit also on the turbulence levels in 
the system. However » a preliminary prediction may be made regarding the 
numerical values of these transfer coefficients through the appropriate 
scaling of water modelling experiments, which were reported in Section 3. 

Upon considering that the conditions of turbulence were comparable 
In the water model system and the pilot scale experiments (indeed the 
water models were designed with this objective in mind) the mass transfer 
coefficients may be scales with the aid of' the following formula: 

N.. = 0.388 (U Tu)N„ 100<RO<2000 Tu^.l3 / (13) 

Nu Re Pr • .* / 

The actual numerical values that one may obtain for the mass transfer 
coefficient ranged from about 0.002-0.05 cm/s which ccmpiire reasonably 
well with those found in the pilot scale, hot metal studies. 

6. DISCUSSION 

In the paper a brief review was presented of recent wor)c concerned 
with the mathematical and physical modelling of transport phenomena in 
gas bubble driven circulation systems of interest in ladle metallurgy. This 
modelling work was augmented by some plant scale measurements. 

The present status of this effort may be summarized as follows: 

(1) Through the statement of the two diir.cnsional turbulent Navier- 
. Stokes equations and the proper boundaries conditions, the use of the 

k-e model provides an excellent representation of the overall flow field 
and gives fairly good guidance regarding the spatial distribution of the 
turbulence parameters. 

(2) It must be realized that the k-e model itself cannot predict the 
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IocaI heat or masf; transfer rates, but that these quantities have to be 
introduced in the form of the l>uundary conditions. In tho absence of 
other reliable information these lioundary conditions have to l>e deduced 
from semi-empirical correlations, which however do contain the computed 
velocity terms, some distance from the interface. It follows that the 
calculation of the local heat or mass transfer rates is neither empirical 
nor entirely fundamentally based, but is some compiromise between these 
two extremes. 

(3) Kxporimenlal measurements usinq water models and Laser Velocimetry 
have been very helpful in providinq a riqorous test of these models. Indeed 
tho above cited comparisons were based on these measurements. On combining 
the known velocity fields with tho direct measurement of the local heat 
transfer rates, through the study of tlio rate at which ice rods melt in an 
agitated water pool valuable insight has been obtained. This work has 
clearly shown that in addition to the local velocities, tho local turbulence 
levels also play a very imi'<ortnnt role in determining the local heat (or mass) 
transfer rates. 

(4) At present only a limited number of plant scale measurements may 
be cited, which may be used for testing tho predictions based on the model. 
However, tho data available indicate that tho combined approach using both 
physical and mathematical models appears to ho very powerful. 

(5) Clearly a number of problems remain unsolved and a great deal of 
furtijcr work \wuld bo desirable in a number of areas, in particular: 

(i) It would bo desirable to obtain a better um'-srstanding of the 
factors that govern tho local h*'at and mass transfer rates between the hath 
and the wall and the bath and a slag cover. 
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(ii) It would be helpful to apply the known spatial distribution 
of the turbulent kinetic energy to develop more comprehensive models 
for deoxidation kinetics. 

(ili) And finally, while the k-c model has proven to be quite 
satisfactory, the development of ah improved model, which has a sound 
theoretical basis should be a future long term goal. 
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List of Symbols 


A 

cross section of area of two phase region 

2 

m 

q 

acceleration due to gravity 

m/s 

11 

right of the liquid column 

m , 

k 

turbulent kinetic energy 

(m/s) , 

p 

pressure 

kg/m. s‘ 

p' 

modified pressure, P-pgz 

kq/ms^ 

m’ls 

Q 

volumetric flow rate 

r 

radial coordinate 

m 

r 

c 

radius of two phase region 

m 

R 

radius of the bath 

m 

U 

velocity 

m/s 

U* 

RMS velocity 


U 

CO 

rising velocity of a single bubble 

m/s 

Z 

axial coordinate 

m 

a 

average void fraction of the gas in 
the two phase region 

2 3 

e 

turbulent energy dissipation 

m /s 

u 

viscosity 

kq/m . s 

V 

kinematic viscosity 

m*/s 

f) 

detisity 

kg/m* 

X 

jet cone angle 



Subscripts 


eff effective 
g gas 
I liquid 
O orifice 
r radial direction 
t turbulent 

z axial direction 0 tangential direction 
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Fig. 1 Relationship between time required for complete 
mixinq and and rate of dissipation of energy 
density i in^^var ious steel processing operations. 
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A comi'iir ison of the expori men tally moasurod and the 
comuutcd velocity profiles for a linear ons velocity 
of i.fi m/s at the orifice, '.b) experimental measurements 
(b) tijooretical predictions. 
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z/H =0.98 


Thu r.vUal dcpomlcnco of thr nbRolutc valxu* of the 
velocity vector - o comparison of the preJictions with 
mocTsurx'mxMils at various axial |>osilionn. Exper im«?ntal 
measxjromcnts for an orifice vilocity xif K>2 cm/s (A) and 
320 rm/s (A). Theoretical predictions for an orifice 
velocity of I f»? cm/s (- . -) and 320 cm/s ( ). 






Fig. 7 The predicted radial variation of the normal velocity 
in- the jet cone region for various axial positions: 

z/H=0.85 cm( ), z=0.60 cm( ), and 

z=0.33 cm ( . ). 
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The CKperimontally determined, normalized radial profile of the 
rms velocity comi>onents. (a) linear orifice velocity 1.62 n/s; 
(b) linear orifice velocity 3.20 m/s; 


Notation; 


for z/ll=U.Bl (V,V) 


Mi). .56(V f). 0.40(0, t),0.30(A A), 0.10(+,#) 









kxIo'^CJ/kg) 
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Fig. 10 Tho radial variation of the turbulent kinetic energy - a 

comparison of tho measurements with predictions at various 
axial positions. Experimental measurements for an orifice 
velocity of 162 cm/s (A and 320 cm/s (A). Experimental 

measurements for an orifice velocity of 162 cm/s( ) 

and 320 cm/s ( ) 
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Fiq. 13 Ice rod diamotor vs. time at z/H=0.3 


f«iTO waiff rsiTH 



















Melting Rate x 10^ (m/s) 


PAGE IS 





z/H 


Fiq. 17 Vortical distribution of the initial melting rate 
at r/R“0.9. 
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